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Abstract 

X 

5j I We present an extensive analysis of rare K decays, in particular of the two neutrino modes 

K^ -^ TT^uV and Kl -^ ir^uV, in the Minimal Supersymmetric extension of the Standard 
Model. We analyse the expectations for the branching ratios of these modes, both within 
the restrictive framework of the minimal flavour violation hypothesis and within a more 
general framework with new sources of flavour-symmetry breaking. In both scenarios, 
the information that can be extracted from precise measurements of the two neutrino 
modes turn out to be very useful in restricting the parameter space of the model, even 
after taking into account the possible information on the mass spectrum derived from 
high-energy colliders, and the constraints from 5-physics experiments. In the presence of 
new sources of flavour-symmetry breaking, additional significant constraints on the model 
can be derived also from the two K^ — > vr^f+f" modes. 



1 Introduction 

As widely discussed in the literature, rare decays dominated by one-loop electroweak 
dynamics offer a very powerful tool to investigate the flavour structure of physics beyond 
the Standard Model (SM) [1]. Among them, the processes K^ -^ vr+z/z/ and K^ —^ tc'^uV 
are certainly a privileged observatory because of the high-level of accuracy achieved in 
their theoretical description [2-7]: future precise measurements of these modes will have 
a non-trivial impact on physics well above the electroweak scale [8]. 

In the last few years precise measurements of flavour-changing neutral-current (FCNC) 
processes in the B sector have severely restricted the parameter space of new-physics 
models, especially in the flavour sector. Moreover, a direct exploration of the physics 
in the TeV range is expected soon with the start of the LHC program. Within this 
context, it is worth to understand if, and at which level, the indirect information which 
could be extracted from rare K decays is still useful. The main purpose of the present 
paper is an attempt to answer this question, within the specific framework of the Minimal 
Supersymmetric extension of the SM (MSSM). 

Several analyses of rare K decays within the MSSM have already been presented in the 
literature, both in the general framework of arbitrary new sources of flavour mixing [9-14] , 
and also in well- motivated scenarios with more restrictive hypotheses [15]. The main 
purpose of all these works has been the identification of the maximal deviations from the 
SM of the two K -^ -kvV rates. Our analysis has a different goal: understanding how 
precise measurements of these observables can be used to discriminate among different 
versions of the MSSM. We will analyse in particular two general frameworks: 

I. The most general version of the MSSM compatible with the Minimal Flavour Vio- 
lation (MFV) hypothesis, as defined in Ref. [16]. 

II. The MSSM with generic new sources of flavour-symmetry breaking, in particular 
with sizable non-MFV trilinear soft-breaking terms in the up sector (A^/), with 
-R-parity conservation and moderate values of tan/? (tan/3 < 30). 

As we will show, in both these frameworks precise measurements of the two B{K — > t^vV) 
are very useful to determine the (flavour) structure of the model. This statement remains 
true even taking into account possible future constraints on the MSSM mass spectrum 
obtained at the LHC, and the refinement of the flavour constraints expected from B 
factories. Within the scenario II, we will show in particular that present constraints still 
allow a large freedom concerning the flavour structure of the Ajy terms. In the presence 
of sizable deviations from the MFV hypothesis in this sector, a key role is played also by 
the two Kl — > 7r°£+£~ modes. 

The paper is organized as follows: in Section |21 we recall some basic formulae for the 
evaluation of rare K decay branching ratios. In Section El we: i) introduce the MFV 
scenario; ii) analyse the expectations of the two B{K — > -kuV) in this framework; iii) 
discuss the consequences of these findings and compare them with the previous literature. 
Similarly, in Section E] we introduce and analyse the consequences for rare K decays of 
the scenario II. The main results are summarized in the Conclusions. 
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2 Basic formulae for rare K decays 

Within the class of models considered here, the supersymnietric contributions to K ^> -nvv 
decays can be described to a good accuracy in terms of a single complex function^ 

w = \Y.w,, (1) 

where Wi are the Wilson coefficients of the following effective Hamiltonian: 



v/2 27r sin^ Qw 



</"" = 5l;£5^ E "'■' --I'd - -^'^^ -^n^i - 75)^, + hx. (2) 



with A = \yus\ = 0.225 ± 0.001. In terms of this function, the two K -^ thju branching 
ratios can be written as 

B{K+^7r^uu) = «:+|1^sm + W^susy|' (3) 

B{KL^7r%u) = KL[lm{WsM + WsvsY)f (4) 

where^ k+ = (5.26 ± 0.06) x lO^^^ and kl = (2.29 ± 0.03) x 10"^°. Defining further 
\q = V*gVqd (where Vij denotes the generic element of the Cabibbo-Kobayashi-Maskawa 
matrix), the SM contribution to the W function reads 

^^. ReAc„ , ReAi .ImAj 

WsM = -^Pu,c + -^^t + ^^^^t ' (5) 

with Xt = 1.464 ± 0.041 [3] and Pu,c = 0.41 ± 0.04 [2,4]. For the numerical values of the 
other input parameters we refer to Refs. [18,19]. 

As discussed in [9-13], among the additional non-standard contributions to the W 
function appearing in the MSSM (W^susy), only those associated with chargino up-squark 
loops and charged-Higgs top-quark loops can compete in size with W^m- A complete 
listing of these contributions can be found in Ref. [13]. 

The effective Hamiltonian necessary to describe the two Kl -^ 7t^£~^£~ modes can in 
general be written as 

^iA5i=i ^ ^^ [ E<(/^)^^^('") + ^^i^^)Q'i^^)] + h-c. , (6) 



^ There are two notable corners of the MSSM parameter space (which will not be analysed in this 
work) where this approximation is not valid: the large-tan /3 scenario with non-MFV right-right mixing 
terms, discussed recently in [14], and models with large violations of lepton-flavour universality, discussed 
in [17]. 

^ The numerical values of the Ki -which encode isospin-breaking and SU{3) violations in the K ^ tt 
matrix elements [7]- have been updated with respect to the previous literature taking into account the 
latest results on the K ^ tt form factors from Ref. [18]. The same comment applies to the Cf in Eq. I^. 



where the hst of potentially relevant operators includes four-quark operators, photon- and 
gluon-dipole operators, and 

Q'jy = S7^(l - 75)^^7/ , Qta = S7^(l - I5)dh,l5^ ■ (7) 

Both within the SM and in the class of supersymmetric models we are considering, the 
direct-CP-violating transition K2 -^ n^i^i^ turns out to be dominated by the contribu- 
tions of Q^ivA- ^^ ^^^^ limit, the corresponding Kl branching ratios can be written as [20] 

B{Kl - 7r°£+r ) = (ci, + c4, + C',^ + C^pc) X 10-l^ (8) 

with the following set of coefficients: 
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Cf„t = -(7.47 ± 0.20) X \as\ x Im w^^y , Cf^i, = (13.9 ± 0.5) x \as 

a;^, = -(1.77 ± 0.04) X \as\ x Im w^y , C^, = (3.27 ± 0.2) x \as\^ , 

C^i, = (2.02 ± 0.10) X [(Im <^)2 + (Im w'.^f] , Qpc ~ , 

C^i, = (1.12 ± 0.05) X [0.43 X (Im wi^yf + (Im wi^^f] , C^p^ = (5-2 ± 1.6) , 

(9) 
Here \as\ = 1.2±0.2 denotes the non-perturbative low-energy constant extracted from the 
experimental value of B{Ks ~^ 7r^i~^i~) [21], and the Wilson coefficients are renormalized 
at the scale /xir ~ 1 GeV. Following the analyses of [20,22], we have assumed a positive 
interference between the long-distance amplitude and the SM short-distance contribution. 
Using the notation of [24], the latter can be expressed as 

. /- \ 27rImAi , „ ,, 27rImAi , , 

{lmw',y)^^ = -yrv x — = - (0.73 ± 0.04) x — ^ , (10) 

(l--r.)sM = -yr- X ^ = (0-68 ± 0.03) x ^ . (11) 

The complete analytic expressions of the chargino up-squark loop contributions to w^yj^ 
can be found in [23]. 

3 The MFV framework 

3.1 Definition of the model 

The MSSM with R parity conservation and generic supersymmetry soft-breaking terms 
has a huge number of free parameters. One of the virtues of B{K -^ ttuu) is that these 
observables are sensitive only to a limited subset of such parameters: those appearing in 
the chargino and up-squark mass matrices, and those which determine the charged-Higgs 
mass. To define the structure of the MSSM we are considering, it is therefore sufficient 
to specify the value of these mass terms. 

The chargino mass matrix in the basis of electroweak eigenstates (wino and higgsino) 

is 

M2 V2Mwsmp 



''^-^v2m:c..p t '• '^'^ 



where the first index of both rows and columns refers to the wino state. Here /i denotes 
the supersymmetric Higgs bilinear coupling, M2 the soft supersymmetry-breaking wino 
mass and tan f3 = Vu/vd the ratio of the two Higgs vacuum expectation values. Note that 
M2 can always be chosen real, without loss of generality. 

Two constraints on the free parameters of M^ could be obtained by measuring chargino 
masses, while independent information could be extracted from the cross sections of vari- 
ous electroweak processes. It is therefore reasonable to assume that the complete structure 
of M^ will be determined, to a good extent, by high-energy experiments.^ The situation 
is very different in the squark sector, where the large number of free parameters does not 
allow a model-independent extraction in terms of high-energy data only. 

The soft-breaking terms appearing in the squark sector are the 3x3 matrices Mg, 
M^, M|, (bilinear terms), and A^, A/) (trilinear terms). Performing in the squark sector 
the same unitary rotations which allow to diagonalize quark mass matrices, 

^L,R ^ ylRqL,R l^CKM = V^Vf (13) 

(i.e. adopting the so-called super-CKM basis), the 6x6 squark mass matrices assume the 
form 



'-uLR ^^^uRR 



clLR dRR 



where 
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mi^^ = V^M^V^Um^ + l-(4M^-Mi)cos2/3, (15) 

M|^^ = V^M2,V^Um^ + l^M|cos2/5sin2^^, (16) 

MIlr = Vu\lAu\l^-cotPfi*m^, (17) 

with m„ = diag(m„, rric, Trit), and similarly for the down sector. For later convenience, we 
also define the couplings 

which parametrize the amount of flavour-symmetry breaking in the up-squark sector (in 
a generic squark basis). 

The flavour structure both of the SM and the MSSM is characterized by a global 
SU{3)^ flavour symmetry, broken only by mass terms and Yukawa interactions [16]. 
Within the SM, the Yukawa interaction is the only source of SU{3)^ breaking. Within 
the MSSM, there are in general several new flavour-symmetry breaking sources, encoded 

^ An additional significant information is obtained by flavour-conserving low-energy experiments, in 
particular by the electric-dipolc-momcnts (e.d.m.) of quarks and leptons, which already provide stringent 
constraints on the possible CP violating phase of fi [25]. 



in the soft-breaking terms. In the following we shall concentrate on the MFV scenario, 
which can be considered as the most restrictive assumption about the structure of these 
additional flavour-symmetry breaking terms. 

According to the MFV hypothesis, the SM Yukawa matrices are the only source of 
breaking of the SU{3)^ flavour symmetry also beyond the SM [16]. Neglecting suppressed 
terms (proportional to high powers of off- diagonal CKM terms and/or light quark masses), 
this symmetry principle implies the following structure for the soft-breaking terms of the 
quark sector (in a generic basis for the electroweak eigenstates) [16]: 



M^ = m^ (^^^ 1 + hYlYu + 62 Y^Yd + k (Y^Y^Yt Y, + Yt Y„ Y^Y, 
M^ = m'l^a^l + hYuYlj , A^; = AY„ (^4! + feeY^Yrf) 

m2(a3l + 65Y,Yj) , A^ = AY, (05! + 67YIY,) , (19) 



M^ 



where a, and 6, are free 0{1) parameters. Redefining these parameters absorbing the 
dimensional factors m and A, this implies for the up squark mass matrix {V = Vckm) 

+m2 + li(4M2,-M|)cos2/5, 
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V: 



w 



'^Ilr = O'^'^u + -^muVmlV'' - cot /?/i*m„ , (20) 



v'd 



in the super-CKM basis. Note that in principle, 04 and b^ can be complex. 

3.2 Scanning of the parameter space 

Employing the parametrizations in Eqs. (j^UI) and (fT^ we have performed a systematic 
scan of the free parameters of the model and analysed the consequences for the two 
B{K — > irvv). Within the MFV framework, the deviations with respect to the SM are 
experimentally undetectable in a good fraction of the parameter space. For this reason, 
we concentrated in particular to identify under which conditions (within this restricted 
scenario) it is possible to generate sizable (detectable) enhancements with respect to the 
SM in the two B{K -^ vrz/z/). The free parameters have been varied in a wide range, 
checking the consistency with tree-level vacuum stability bounds [26] , direct experimental 
constraints on squark and gaugino masses (see Table P) and existing constraints from 
precision measurements (both in the electroweak and in the flavour sector). The scan has 
been performed using the adaptive integration routine VEGAS [27], to search for maximal 
effects and to deal most effectively with invalid regions in parameter space. 



Ranges (GeV) 


Exp. bounds (GeV) 


ai_3 


[0, 1000] 


M+ > 94 


M,,,,,,j > 250 


04 


0, 3000] 


M^o (M^o) > 46 (63) 


M4>96 


l/^l 


[0, 500] 


M^o (M^o) > 100 (116) 


Ml >89 


M2 


0, 3000] 







Table 1: Ranges adopted in the scan for the free parameters of the MFV framework 
relevant for K — ^ iri'T'. Experimental lower bounds on squark and gaugino masses [28]. 



The chosen ranges of the relevant parameters are reported in Table ^ The parameter 
05 and all the bi have a negligible effect on the considered branching ratios and can be set 
to zero without loss of generality. The MFV implementation of vacuum stability bounds 
then assumes the simple form |a4p < 3(a^ + a"^)- Also a^ has a very small impact, but 
has to be taken sufficiently large to generate down squark masses large enough to pass 
the experimental bounds. A similar comment applies to the gaugino mass parameter Mi 
in connection with the experimental bounds on neutralino masses. For definiteness, we 
have taken Mi = 500 GeV. As far as tan/3 is concerned, we have fixed it to the reference 
value tan/3 = 2: as we will discuss in the following, larger/smaller values of tan/3 lead 
to smaller /larger effects in B{K -^ ttuu).^ The impact of the phases of /i and 04 has 
also been explicitly studied, and found to be quite small (we will come back to this point 
later), so we considered only real parameters in the final scan. In any case, the different 
sign combinations of fi and 04 allow to keep track of these phase effects. Finally, in 
order to disentangle the various contributions, we have set Mh+ > 1 TeV, such that the 
charged-Higgs top-quark loops decouple. As chargino box effects can be safely neglected 
for a lightest slepton mass above 300 GeV, in practice we are left only with the chargino 
up-squark Z penguin. 

The main results of the scan are illustrated in Figures C] and |21 As shown in Figured 
the maximal enhancements of the branching ratios are closely correlated to the minimal 
values for the lightest stop and chargino masses.^ This correlation is very useful to de- 
termine the flavour structure of the model. For instance, if B{K^ -^ vr+z/z/) was found to 
be more than 10% above the SM expectation and the lightest stop and chargino masses 
were both found to be above 130 GeV, with a charged Higgs mass above 1 TeV, then one 
could exclude the MFV scenario. As expected, the MFV hypothesis predicts also a strict 
correlation between B{K~^ -^ tt^uu) and B{Ki^ — >■ ir^uu). We flnd in particular 



R{K+ -^ TT+z/z/) = (0.965 ± 0.008) x R{Ki^ -^ n^i^iy) 



(21) 



^ Low values of tan/3 are strongly constrained by the experimental constraints on the Higgs sector, 
and in particular by the lower bounds on the lightest neutral Higgs [29] . We have explicitly checked that 
for tan/3 = 2 (and even slightly below) we are compatible with these experimental bounds. This happens 
mainly because we allow sizable trilinear soft-breaking terms in the up sector. 

^ The variation of density in these plots is not strictly correlated to the density of the underlying 
parameter space: by construction, we scan with more points the regions with larger effects. 



E 1. 




Figure 1: Regions in the rrif - m^ plane (lightest stop and chargino masses) allowing en- 
hancements oi B{K~^ -^ T^^uu) of more than 11% (yellow/light gray), 8.5% (red/medium 
gray) and 6% (blue/dark gray) in the MFV scenario, for tan/3 = 2 and Mh+ > 1 TeV [the 
corresponding enhancements for B{Ki^ -^ tt^uu) are 15%, 12.5% and 10%, respectively, 
see Eq. (ED)]. 



R{K ^ /) = BiK ^ f)/BiK ^ /)sM , (22) 

in the region of maximal enhancements (i.e. 10% to 16% for the neutral mode). In 
principle, the relation (PT|) would allow the best test of the MFV hypothesis. However, the 
experimental challenges of the i^L -^ ir'^vv mode make the correlation between B{K~^ — >■ 



7r"^z/z/ , m 



and m^ outlined in Figure Q a more useful test for the near future. 



In Figure |21 we present a more detailed analysis of the parameter-space region with 
enhanced B{K — >■ vrz/P), showing the two-dimensional projections on the most significant 
planes. In this case is even more evident the key role of a precision measurement of B{K -^ 
vrz/z/) in selecting a well-defined region of the model, or in constraining its structure. As 
can be noted, an important role is played by the parameter a^. sizable enhancements of 
B{K — >■ T^uu) can occur only for large enough values of this parameter. The reason of this 
effect can be traced back to the enhancement mechanism discussed in Ref. [11]. Indeed, 
even within the MFV scenario one generates non- vanishing left-right flavour-mixing terms 
in the squark basis of Refs. [10,11]. In particular, the double mass-insertion combination 
which controls possible enhancements in B{K —>■ ttuu) [11,13] assumes the form^ 



(23) 
We denote by {SYnJij the flavour-mixing couplings of Eq. p8|l in the squark basis of Refs. [10, 11]. 
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Figure 2: Correlations of the most significant MSSM-MFV parameters corresponding 
to different enhancements of E{K^ -^ n^uu) (identified by the color/gray-scale as in 
Figure HJ . 



and thus grows with a^. The expression ()23p also shows that: i) a possible CP-violating 
phase of a^ has a negligible impact if fi is approximately real (as implied by the stringent 
e.d.m. bounds); ii) maximal effects are obtained for /x and 04 of opposite signs; iii) small 
values oita.n/3 enhance the magnitude of (5rl)32(^Kl)3i ^^^^^ further. We finally note that 
the contribution of fi is always subleading with respect to the one of 04 (and thus the 
sensitivity to tan/5 is quite mild) since large values of /i lead to a strong suppression of 
the loop function (because of heavy higgsino masses). 

As far as the constraints from other observables are concerned, the MFV hypothesis 
automatically implies small non-standard effects in FCNC observables such as ex and B -^ 
Xg^. Moreover, these observables are sensitive to different combinations of free parameters 
with respect to the two B{K — > ttuu). As a result, the requirement of consistency with 
existing FCNC data does not have a perceptible impact on Figures [T] and El As pointed 
out in Ref. [15], an important constraint is obtained from flavour-conserving electroweak 
observables and, in particular, from Ap. The correlation between Ap and B{Ki^ — »• tt^vv) 
is illustrated in Figure El The effect for B{K^ —^ tx^vv) is almost identical, provided one 
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Figure 3: Correlation between Ap and B{Ki^ 
tan/? = 2 and Mh+ > 1 TeV). 



TT 



'^uu) in the MFV framework (with 



scales the enhancement according to Eq. (PT|) . 

3.3 Comparison with previous literature 

We conclude the MFV analysis with a comparison of our results with those obtained in 
a similar framework in Refs. [15] and [30]. First of all, we recall that for simplicity the 
numerical results in Figures Q and El have been obtained in the limit where the charged- 
Higgs contribution can be neglected. This approximation can easily be removed: charged- 
Higgs top-quark loops always induce a constructive contribution to the W function, which 
depends only on tan/? and Mh+ (see e.g. Ref. [11]). For tan/? = 2, the maximal effect 
amounts to an additional ~ 5% enhancement of the two B{K — > vrz/z/) (whose limiting 
factor is the Mh+ ^ 300 GeV bound derived from B -^ X^j). As a result, we conclude 
that in the MFV framework 



R{K+ -^ TT+z/z/) < 1.20 , R{Kl -^ 7r%u) < 1.25 



(24) 



This conclusion is in contradiction with respect to the claim of Ref. [15] that no sizable 
enhancement of the two B{K — > vrz/z/) is possible in the MFV framework. The origin of 
this difference is twofold: 

• We have considered a more general definition of the MFV framework (the most gen- 
eral definition compatible with the renormalization group). Having a more general 
structure for the soft-breaking terms, within our scheme the effects in B{K -^ vrz/z/) 
are less severely constrained by the existing constraints on other observables. 

• We have adopted a different strategy concerning the fit of the CKM matrix. In our 
analysis we have implicitly assumed that within the MFV framework the determi- 
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nation of the CKM matrix is not affected by the presence of new physics. Thanks 
to the recent precise results of B factories, we now know that this assumption is 
an excellent approximation (see e.g. [31]). At the time of Ref. [15], the available 
experimental information was less precise and the authors decided to perform a 
non-standard CKM fit using observables sensitive to SUSY corrections. This fact 
introduced spurious correlations between the genuine SUSY effects and indirect ef- 
fects associated with the CKM determination. In particular, the predictions of the 
two B{K —>■ nuu) turned out to be suppressed because of smaller effective values of 
\t, and not because of suppressed loop functions. 

In conclusion, we find that in the MSSM realization of the MFV hypothesis the two 
B{K — > irui?) can saturate the model-independent bounds of Ref. [30]. This implies that 
a measurement of one of the two B{K — >■ irui?) consistent with Eq. P^ does not allow 
to distinguish a generic MFV model from the MSSM. On the other hand, as illustrated 
by Figure ^ stringent tests of the model can be performed combining B{K -^ txvv) and 
sparticle mass measurements. 

4 The general framework 

As anticipated, within the MSSM there are in principle several new (non- Yukawa) sources 
of flavour-symmetry breaking. Some of them are highly constrained by precise data on 
various rare processes. However, we are still far from being able to conclude that the 
MFV hypothesis, namely the absence of new sources of flavour-symmetry breaking, is the 
only viable option. 

Focusing the attention on the quark sector, we can identify five independent sources of 
flavour-symmetry breaking in the matrices Mg, M^, M|,, Ajy, and A^). Their non-trivial 
parts (terms not proportional to the identity matrix) introduce breaking terms of the 

GT::L = SUi3)Q, X SUi3)u, X SU{3)n, (25) 

subgroup of SU{3)^ transforming respectively as 

M|~(8,l,l), M2,~ (1,1,8), As ~ (3,1,3), (26) 

Af;~(3,3,l), M2,~ (1,8,1). (27) 

In the literature, there exist several phenomenological analyses of such terms, typi- 
cally expressed as upper bounds on the mass- insertion couplings ((5^|g )jj, defined as in 
Eq. (fTH|) . In particular, ex and AM^ imply stringent bounds on all the 1-2 down-type 
mass insertions (limits in the 10~^-10~^ range for squark masses below 500 GeV) [32]; 
bounds in the 10~^ range for all the 1-3 down-type couplings follow from AM^^ and 
Acp{Bd — >■ J/ipK) [33]; bounds in the 10^^-10^^ range are derived on the 2-3 down-type 
couplings from B -^ Xs'y [34] (which strongly constrains LR terms) and recently also by 
AMb^ [34] (which is also sensitive to LL and RR terms). All these stringent phenomeno- 
logical limits have been derived analysing the impact of gluino-mediated amplitudes. As a 
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Charginos: 


/i = 500 ± 10 GeV 


M2 = 


= 300 ±10 GeV tan/? = 2-4 


Up-squarks: 


Muj, = 600 ± 20 GeV 


M,,, 


= 800 ± 20 GeV A^ = 1 TeV 


Other mass terms: 


Ml = 500 GeV 




= M; = M3 = Mh+ = 2 TeV 



Table 2: Basic choice of the flavour-conserving parameters used in Figures EHHl 

consequence, the constraints concern only down-type mass matrices, or the three flavour- 
symmetry breaking structures in Eq. (j^ . The bounds on the up-type soft-breaking 
terms, derived from chargino amplitudes, are substantially weaker (see e.g. Ref. [35] for a 
recent analysis in the 2-3 sector). 

Interestingly, within the SM the only large breaking of the flavour symmetry appears in 
the up-sector, or in the up-type Yukawa coupling transforming as (3, 3, 1). It is therefore 
quite natural to conceive supersymmetric scenarios where the three flavour-symmetry 
breaking structures in Eq. ()2fi|l are very small (in agreement with observations) and sizable 
non-minimal breaking terms appear only in the up- type structures in Eq. ()27|1 . especially 
in the (3, 3, 1) sector. As we shall show in the rest of this section, this scenario is perfectly 
compatible with all existing constraints on B and K physics, and rare K decays are the 
most useful tools to probe it in the future. 

4.1 Definition of the model 

We consider a non- minimal scenario where Mg, M|,, and M^ have an approximate MFV 
structure, while A(/ contains sizable non-minimal flavour-breaking terms of 0{X), where 
A ~ 0.22 is the Cabibbo angle. More precisely, we assume that in the super-CKM basis 
M|^^ and M?j:j^ have the form in Eq. (^0]) . while 

MLii=[A-cot/3/i*]m„, (28) 

where 

(A)33 = Ao, |(A),^,|<AAo. (29) 

This non-minimal structure is naturally consistent with vacuum stability bounds [26]. We 
note that structure of the type (j^Hj) also for Mj (or Ad), with lUu ^->- nirf, is generally 
consistent with the existing bounds from gluino-mediated FCNC amplitudes (given the 
smallness of m^). However, in order to isolate the effects induced by A^/, in the following 
we concentrate on the case where A^ is aligned to the corresponding Yukawa matrix. 

In order to compare the sensitivity of various FCNC observables to the new sources 
of flavour-symmetry breaking, we fix the flavour- conserving parameters of the model to 
some reference values, and study the dependence of the observables on the Aij terms. 
This way we simulate somehow a post-LHC scenario, where most of the flavour- conserving 
parameters of the model are known, because of the progress at the high-energy frontier, 
while precision measurements of rare decays can be used to determine the flavour structure 
of the model. 
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4.2 Numerical analysis 

The basic choice of the flavour- conserving parameters is shown in Table 121 For simplic- 
ity, a high scale (2 TeV) has been chosen for all the parameters which play a minor role 
in chargino up-squark amplitudes. The results depend very little from this choice, pro- 
vided the minimum value of down-squark/slepton masses is above 500 GeV. We have 
also assigned a small error to the key mass parameters of the chargino up-squark/sector, 
in order to investigate the sensitivity of flavour- changing observables to the precision on 
mass terms. Given the smallness of the Yukawa couplings of the first 2 generations, the 
only relevant Aij terms are ^13, A23, and ^33 = Aq. As in the previous section, we have 
varied the free parameters of the model checking the consistency with direct experimental 
constraints on squark and gaugino masses and the existing constraints from rare processes 
(in particular the B — >■ Xg'y one). 

The numerical results thus obtained for the two K —>■ ttuu modes, in comparison with 
Bd -^ Xsi'^i~, Bd -^ /i^yU~, AMb^, and the CP- violating observables ex and Acp{Bd — >■ 
J/ijjK), are shown in Figures EHHl The main features resulting from this numerical study 
can be summarized as follows: 

• The non-standard effects induced by these chargino-mediated amplitudes, in the 
presence of non-MFV up-type A terms, are maximal in the two K —>■ ttuu decays. 
The dominance of i^ ^ nui? holds in comparison with other K- and i?-physics 
FCNC amplitudes, both in CP-conserving (Figure llj and in CP-violating observ- 
ables (Figure E}. 

Note that the non-standard effect in B{K —>■ ttuu) is not necessarily an enhancement 
with respect to the SM. For instance, for A13 ^ one has B{K^ -^ -k^vv) < 
B{K^ -^ 7r+z/z/)sM, while in the same region of parameter space AM^^ receives a 
~ 4% positive correction (see Figure |31). However, the important feature emerging 
from our analysis is that B{K^ — > ir'^ui?) has a much stronger dependence on Ai^ 
with respect to AMb^ (and the other i?-physics observables). 

• As can be noted from Figures EHSl despite the fact that the flavour conserving 
parameters of the model are completely fixed, the predictions for the FCNC observ- 
ables are still affected by a large uncertainty. This is due to the free parameters 
hidden in phases and moduli of A13 and A23 (see Figure EI). For this reason, a 
complete determination of the model requires several precision measurements in the 
FCNC sector. Ideally, the two K — >■ nuu rates and two clean FCNC observables in 
the B system (such as the rates of the two Bs^d — ^ /U"*"/!" modes). 

• As can be seen in FigureEl even setting A23 = 0, the dependence of B{K -^ tipt)) on 
Ai3 is quadratic, contrary to the approximate linear dependence of the other FCNC 
observables. This quadratic dependence, which is the main reason of the dominance 
of non-MFV chargino-mediated amplitudes in K -^ ttuu, can be explained in terms 
of the double- mass insertion mechanism of Ref. [11]. 
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Figure 4: Dependence of various FCNC observables (normalized to their SM value) on the 
up-type trilinear terms A13 and A23, varied according to Eq. ()29|) . The flavour- conserving 
parameters of the model are flxed as specifled in Table El Upper plot: B{K~^ — >■ tt'^ui') 
(blue/dark gray), B{B(i — > yU^/U^) (red/gray lower- region) , AMb^ (green/gray upper- 
region) as a function of A13. Lower plot: same observables as in the upper plot, with the 
superposition of B{Bd -^ Xgi^i") (light green/light gray) plotted as a function of A23 
(instead of A13). 
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Figure 5: Anatomy of the SUSY contributions to B{K~^ — > vr+z/z/), B{B(i -^ yW"*"/^") and 
AMb^ changing the range of variability of A13 and ^423. Notations and other conventions 
as in Figure El 



To better understand this effect, we denote by Aij (Aij) the A terms in the super- 
CKM basis (the squark basis of Refs. [10, 11]). It is easy to realize that if A23 = 
and A13 7^ in the super-CKM basis, then both A13 and A23 are not vanishing 
because of the non-diagonal entries of the CKM matrix. In particular. 



A. 



23 



V2jA,3 



V2iAr3 



(30) 



A similar mechanism occurs if ^413 = and A23 7^ 0. Thus, barring fine-tuned 
scenarios, there is always a sizable effective double mass-insertion coupling, given 

by 

{Srl)12 i^RL)31 ^ ^23Al, ^ Vu\A23\' + 1^2l|Al3p + ^23^t3 (31) 

Since the double mass-insertion enhancement is effective only in the kaon system 
(because of its double-CKM suppression [11]), rare K decays are naturally the best 
probe of any non-MFV structure in Ajy. 

A further check of the dominance of the double mass insertion mechanism is shown 
in the lower plot of Figure El for the CP- violating channel i^L "^ vr^z/z/. Here both 
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Figure 6: Impact of chargino- mediated amplitudes on CP- violating observables. Up- 
per plot: comparison of ex (bordeaux/dark gray), Aq^[B(1 -^ J/ipK) (violet /gray) and 
B{Ki^ —>■ -K^vv) (light blue/light gray). Lower plot: dependence of the non-standard con- 
tributions to B{Ki^ -^ TT^ui?) on the CP-violating phases of A13 and A23 in the region of 
destructive interference between SUSY and SM contributions (light-blue/light-gray dots 
correspond to no constraints on the phases). The trilinear terms are varied according to 
Eq. (j211); all the other supersymmetric parameters are fixed as in Table El 
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Figure 7: Comparison of the decoupling of non-standard contributions to ex (bor- 
deaux/dark gray) and B{Ki^ — >■ 'k^uv) (light blue/light gray) in the limit of heavy su- 
persymmetric particles. The scatter plots are obtained varying Ais and ^23 according to 
Eq. (j2ni), ^ur in the interval 200-1000 GeV, and fixing all the other mass parameters as 
in Table 121 The horizontal axis denotes the lightest up-type squark mass. 



^23 and Ai^ are non vanishing, and regions corresponding to specific choices of 
their CP- violating phases are outlined with different colors (gray scales). As can 
be noted, to a good approximation the non-standard effect depends only on their 



relative phase and -consistently with Eq. (|M|)- 
imaginary. 



is maximal when A23A]^3 is purely 



At first sight, it is quite surprising that the non-MFV left-right mixing terms in 
Eq. (j2ni) are not excluded by the precise data on B{B — > Xs7) and have a marginal 
impact also on B{B — > Xsi^i~) (see Figure EJ. However, this fact can easily be 
understood by noting that a non- vanishing b ^ s'j amplitude (generated by effective 
operators of the type hRa^ySiF^^" or bicr^^SnF^'^) requires: 

i. odd number of chirality flips in the down sector; 
ii. odd total number of chirality flips summing up, down and flavour-independent 



16 



5-lxLO 
4.8xL0 
4JxL0 
4.2XLU 
3-9xL0 
3.6x10 
3-3xL0 
3.0x10 
2.7x1 



BR(K|^^7r ee) vs BR(K|^^m'A') 




p(37.5 ± 7.5)= 
2.5 ± 7.5)= 
7=(67.5 ± 7.5)= 
"H8^"5 ± 7.5)= 
"^(97.5 1 7.5)= 



0.0 



2.0x10" 



4.0x10" 



6.0x10" 



8.0x10" 



1.0x10" 



TT^e^e 



Figure 8: Correlation between B{Ki^ -^ tt^uu) (horizontal axis) and B{Kl 
(vertical axis) in the MSSM scenario with non-minimal ^13 and A23, varied according 
to Eq. ()29|) (other parameters are fixed as in Table |21). Different colors (gray scales) 
correspond to different values of the CKM phase 7. 



terms. 

This implies that the up-type left-right mixing terms in Eq. (j29p can contribute to 
the 6 ^ 57 transition only via amplitudes which have at least three chirality flips 
(the up-type trilinear term, one left-right mixing in the down sector, and a third 
SU{2)l breaking term in order to recover the total helicity- violating structure). 
Since each chirality flip is associated with an insertion of the SM Higgs vev, this 
structure implies a strong suppression. 

A phenomenological check of this statement is provided by the loose bounds on 
(^Kl)32 extracted from B{B -^ Xg'j) (see e.g. Refs. [35,36]). 

An important property of the chargino-mediated non-MFV contributions is their 
slow decoupling in the limit of heavy superpartners in K ^ iruu decays (~ m. 



SUSYy'; 



compared to the fast decoupling in AF = 2 observables (~ ""^susy)- This property, 
which has been discussed in detail in Ref. [12], is clearly illustrated in Figure [7| 
in the case of CP-violating observables (a completely similar scenario holds in the 
CP-conserving case). This property provides another natural explanation of why we 
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can still hope to observe sizable deviations from the SM in rare K decays, despite 
the absence of non-standard effects in AM^^^, y4cp(-Bd — > J/ipK), and ex- 



• A well-defined prediction of this non-MFV scenario is a strict correlation between 
the non-standard contributions to Ki^ —^ n^uu, Kl — ^ TT^e^e~ , and Kl — > vr^/i+yU". 
As shown in Figure the relative size of the effect is larger in the Ki^ -^ tt^uu 
case; however, observable deviations from the SM are expected also in the two 
Kl — >■ 7T^i~^i~ modes. This correlation is a consequence of the Zsd effective vertex 
common to the three decay modes, which encodes the dominant non-SM contribu- 
tions in the limit of heavy sleptons. The detailed structure of the correlation is 
almost independent of the flavour-conserving parameters of the MSSM, for suffi- 
ciently heavy sleptons, but has a significant dependence of the CKM phase 7 (as 
shown in Figure |H1). 

If at least two of the rare Kl channels could be measured with good precision, their 
correlation within this non-MFV scenario would provide a striking signature of the 
model. 

5 Conclusions 

The determination of the fiavour structure of the MSSM -as well as of any TeV extension 
of the SM- is one of the components of the so-called inverse problem [37], which hopefully 
we will face soon with the start of the LHC program. While the direct exploration of the 
TeV scale will reveal some of the features of the new-physics scenario, the LHC program 
alone is unlikely to completely determine the structure of the new underlying theory. 
This statement is particularly true in the case of the fiavour structure of the theory, 
whose mo del- independent determination requires new high-precision measurements also 
at low energies. 

In this paper we have provided a quantitative illustration of this problem, analysing 
the impact of the two K -^ ttuV modes in determining the fiavour structure of the 
MSSM. We have analysed the expectations for the branching ratios of these modes in 
two representative classes of the MSSM (as far as the fiavour structure is concerned): the 
Minimal Flavour Violation scenario [16] and a scenario with Af/ terms not aligned with 
the corresponding Yukawa coupling. In both these frameworks precise measurements of 
the two B{K -^ TTuV) turn out to be key ingredients to determine the structure of the 
model. 

Within the MFV scenario, the deviations from the SM in the two B{K —>■ ttuV) are 
naturally small (typically within 10%), but could saturate the model- independent bounds 
of Ref . [30] . We have outlined some clean correlations between the possible non-standard 
effects in the two B{K -^ iruV) and the values of stop and chargino masses which provide 
a distinctive signature of this scenario, or an efficient method to test the implementation 
of the MFV hypothesis in the MSSM. 

The situation is certainly more interesting in the case of non-MFV up-type trilinear 
terms. Within this scenario, which is well motivated and perfectly compatible with all 
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existing constraints from B and K physics, the two B{K — > nuV) could receive 0{1) cor- 
rections with respect to the SM. We have indeed demonstrated that these rare K decays, 
together with the two Kl —>■ n^i~^i~ modes, represent the most sensitive probe of any 
misahgnement between Au and the corresponding Yukawa couphng. The precise mea- 
surement of these rare decays is therefore a necessary ingredient for a model-independent 
reconstruction of the flavour structure of the MSSM soft-breaking terms. 

Acknowledgements 

We thank Andrea Brignole and Simone Pacetti for useful discussions. This work is par- 
tially supported by IHP-RTN, EC contract No. HPRN-CT-2002-00311 (EURIDICE). The 
work of CS. is also supported by the Schweizerischer Nationalfonds; S.T. acknowledges 
the support of the DFG grant No. NI 1105/1-1 and DFG within SFB/TR-9 (Computa- 
tional Particle Physics). 

References 

[1 



[2 

[3 

[4] 
[5] 

[6 

[7 
[8 

[9 

[lo; 



See e.g. G. Isidori, Int. J. M od. Phys A17 (2 002) 3078 ,hep-ph/01 10255 1; A nnales 
Henri Poincare 4 (2003) S97 hep-ph/030lT59|; A.J. B uras, |hep-ph/0505175| L. Sil- 



vestrini, 



hep-ph/05 10077} T. Hurth, hep- ph/051 1280 and references therein 



A. J. Buras, M. Gorbahn, U. Haisch and U. Nierste, |hep-ph/0508165| 
|hep-ph/0603"079 

G. Buchalla and A. J. Buras, Nucl. Phys. B548 (199 9) 309 |hep- ph/9901288j ; 
M. Misiak and J. Urban, Phys. Lett. B451 (1999) 161 |hep-ph/9901278p 



G. Isidori, F. Mescia and C. Smith, Nucl. Phys. B718 (2005) 319 ,hep-ph/0503107| . 

A. F. Fal k, A. Lewandowski and A. A. Petrov, Phys. Lett. B505 (2001) 
107 |h ep-ph/0012099|; M. Lu and M. Wise, Phys. Lett. B324 (1994) 461 
|hep-ph/9401204j . 

G. Buchalla and G. Isidori, Phys. Lett. B440 (1998) 170 |hep-ph/980650ll . 

W.J. Marciano and Z. Parsa, Phys. Rev. D53 (1996) Rl. 

D. Bryman, A. J. Buras, G. Isidori and L. Littenberg, Int. J. Mod. Phys. A21 
(2006) 487 hep-ph/050517H; and references therein. 



Y. Nir and M.P. Worah, Phys. Lett. B423 (1998) 319 | hep-ph/9711215| . 

A.J. Buras, A. Romanino and L. Silvestrini, Nucl. Phys. B520 (1998) 3 
|hep-ph/9712398| . 



[11] G. Colangelo and G. Isidori, JHEP 09 (1998) 009 |hep-ph/98 08487| 

19 



[12; 
[is; 

[14 

[is; 

[16 

[17; 
[is; 

[19 

[2o; 

[21 

[22; 
[23; 

[24 

[25; 

[26 
[27; 
[28; 
[29 
[30 
[31 



A.J. Buras et a/., Nucl. Phys. B566 (2000) 3 |hep-ph /990837l] . 

A. J. Buras, T. Ewertli, S. Jager and J. Rosiek, |hep-ph/0408142[ 

G. Isidori and P. Paradisi, |hep-ph/Q601094[ 

A. J. Buras et ai, Nucl. Phys. B592 (2001) 55 |hep-ph/0007313|. 

L. J. Hall and L. Randall, Phys. Rev. Lett. 65 (1990) 2939; G. D'Ambrosio, 
G. F. Giudice, G. Isidori and A. Strumia, Nucl. Phys. B645 (2002) 155. 

Y. Grossman, G. Isidori and H. Murayama, Phys. Lett. B588 (2004) 74. 



E. Blucher et al, hep-ph/05 12031^ 



Ihttp : //www . slac . Stanford . edu/xorg/hf ag[ http : //utf it . romal . inf n . it| 



G. Buchalla, G. D'Ambrosio and G. Isidori, Nucl. Phys. B672 (2003) 387 
|hep-ph/0308008|; G. Isidori, C. Smith and R. Unterdorfer, Eur. Phys. J. C36 

(2004) 57 hep-ph/040412 7|. 

J.R. Batley et al. [NA48/1 Collaboration], Phys. Lett. B576 (2003) 43 
|he p-ex/0309075i . 

S. Fri ot, D. G reynat and E. De Rafael, Phys. Lett. B595 (2004) 301 
ihep-ph/0404136| . 

P. L. Cho, M. Misiak and D. Wyler, Phys. Rev. D54 (1996) 3329 |hep-ph/960l"360| . 

A.J. Buras, M.E. Lautenbacher, M. Misiak and M. Miinz, Nucl. Phys. B423 (1994) 
349 |hep-ph/9402347|. 



S. Pokorski, J. Rosiek and C. A. Savoy, Nucl. Phys. B 570 (2000) 81 
| hep-ph/9906206j . 

J. A. Casas and S. Dimopoulos, Phys. Lett. B387 (1996) 107. 

G.P. Lepage, J. Comput. Phys. 27 (1978) 192. 

S. Eidelman et al. [Particle Data Group Collaboration], Phys. Lett. B592 (2004) 1. 

A. Sopczak [ALEPH Collaboration] , |hep-ph/0602136[ 



C. Bobeth et a/., |hep-ph/0505110| 

F. J. Botella, G. C. Branco, M. Nebot and M. N. Rebelo, Nucl. Phys. B725 (2005) 
155 [hep-ph/0502133; M. Bona et al. [UTfit Collaboration] , _hep-ph/0509219. 



[32] M. Ciuchini et al, JHEP 9810 (1998) 008 |hep-ph/9808328| . 

20 



[33] D. Becirevic et ai, Nucl. Phys. B634 (2002) 105 |hep-ph/0112303 
[34] M. Ciuchini and L. Silvestrini, |hep-ph/0603TT^ 
[35] Z. Xiao, F. Li and W. Zou, |hep-ph/0603120' 



[36] M. Misiak, S. Pokorski and J. Rosiek, Adv. Ser. Direct. High Energy Phys. 15 
(1998) 795 jhep-ph/ 97034421 . 



[37] N. Arkani-Hamed, G. L. Kane, J. Thaler and L. T. Wang, hep-ph/05 12190] 



21 



